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Abstract 

We compute the static-light baryon spectrum by means of Wilson twisted mass lattice QCD us- 
ing A'^^ = 2 flavors of sea quarks. As light ujd valence quarks we consider quarks, which have the 
same mass as the sea quarks with corresponding pion masses in the range 340 MeV < mps ^ 525 MeV, 
as well as partially quenched s quarks, which have a mass around the physical value. We consider 
all possible combinations of two light valence quarks, i.e. A, S, S and baryons corresponding 
to isospin / G {0 , 1/2, 1} and strangeness G {0 , —1 , —2} as well as angular momentum of 
the light degrees of freedom j G {0 , 1} and parity P G {+ , — }. We extrapolate in the light ujd 
and in the heavy h quark mass to the physical point and compare with available experimental 
results. Besides experimentally known positive parity states we are also able to predict a number 
of negative parity states, which have neither been measured in experiments nor previously been 
computed by lattice methods. 



1 Introduction 



In this work we report on a lattice computation of the spectrum of b baryons made from a heavy 
b quark and two hght quarks, which are u, d and/or s. 

Experimentahy five b baryon states have been observed. While has first been detected quite 
some time ago, S;,, S^, Hf, and ^lh have only been discovered recently [U [21 [3l IU [5] . For the mass 
of rjf, there are two different results, which are not in agreement. 

On the theoretical side there are a number of lattice studies of the spectrum of b baryons. Some 
of these consider static heavy quarks [6l O [HI [9l [lOl [11] using Heavy Quark Effective Theory 
(HQET) (cf. e.g. [121 [13]), while others apply heavy quarks of finite mass [l4l[l5l[T6] mainly by 
means of Non-Relativistic QCD (cf. e.g. [E]). For a recent review of lattice results on b baryon 
masses cf. [18]. 

In this work we treat the b quark in leading order of HQET, which is the static limit. In this limit 
there are no interactions involving the spin of the heavy quark, i.e. states are doubly degenerate. 
Therefore, it is common to label static-light baryons by integer spin/angular momentum j and 
parity V of the light degrees of freedom. For the two light quarks we consider all possible 
combinations of u, d and s, i.e. further quantum numbers are strangeness S and isospin I. We 
use Nf = 2 flavors of dynamical quarks and study various ensembles with corresponding pion 
masses down to ~ 340 MeV. Our lattice spacing a ~ 0.079 fm is rather fine and we use the Wilson 
twisted mass formulation of lattice QCD at maximal twist, which guarantees automatically 0{a) 
improved spectral results. We compute all five experimentally known b baryon states. We also 
make predictions for H^, which has not yet been observed, as well as for a number of negative 
parity static-light baryons, which have neither been measured experimentally nor been computed 
by lattice methods. 

The next-to-leading order of HQET, which removes the degeneracy with respect to the heavy 
quark spin, is ©(l/mg), where mg is the mass of the heavy quark. This correction is expected 
to be relatively small for b baryons, e.g. experimentally the mass difference between S;, and 
is only around 21 MeV. Lattice methods to evaluate such l/mg contributions have been estab- 
lished and tested in quenched studies of B mesons [TH [201 [ZD [22] . We intend to explore these 
contributions using lattice techniques subsequently. An alternative way to predict the spectrum 
of b baryons is to interpolate between charmed baryons, where the experimental spectrum is 
rather well known, and the static limit obtained by lattice QCD assuming a dependence as 
l/mg. Thus the splittings among b baryons should approximately be mc/mi) ~ 1/3 of those 
among the corresponding c baryons. 

We try to determine the b baryon spectrum as fully as possible, i.e. we consider all possible 
light flavor combinations corresponding to 5 G {0 , —1, —2} and / G {0, 1} as well as both 
parity V = + and V = —. This will help the construction of phenomenological models (cf. 
e.g. |23j). might contribute to resolve open experimental issues (e.g. the above mentioned mass 
discrepancy for il;,) and also provide valuable input for future experiments. 

This study is in many aspects similar to our recent computation of the static-light meson spec- 
trum [241 [25] . Preliminary results have already been reported in conference proceedings |26j . 

The paper is organized as follows. In section [2] we briefly recapitulate our lattice setup, which 
is discussed in more detail in |24j . In section [3] we discuss static-light baryon trial states. 
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corresponding correlation matrices and how we extract the static-Hght baryon spectrum from 
these matrices as well as our extrapolation procedure to the physical u/d quark mass. In section^] 
we interpolate between our static-light lattice results and experimental results for c baryons, to 
account for the finite mass of the h quark. We conclude with a brief summary and an outlook 
in section O 
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2 Lattice setup 



In this work we use the same setup as for our recent computation of the static-hght meson 
spectrum. For a more detailed presentation we refer to [24:\ 125] , 

We use Nf = 2 flavor gauge field configurations generated by the European Twisted Mass 
Collaboration (ETMC). The gauge action is tree-level Symanzik improved j27j . 

Sg[U] = E E (1) 

with bo = 1 — 8bi and hi = —1/12. The fermionic action is Wilson twisted mass (cf. [28 t [29 t [30 l 

EH), 

Sf[x,X,U] = a^'^x{x)(^D^^ + ifiq-f5T3^x{x), (2) 

X 

where 

= ^(7m(v^ + V;) -aV;^V^) +mo, (3) 

and V* are the gauge covariant forward and backward derivatives, itiq and ^Uq are the bare 
untwisted and twisted quark masses respectively, rs is the third Pauli matrix acting in flavor 
space and x = {x^^\x^'^^) represents the quark fields in the so-called twisted basis. The twist 
angle uj is given by tan(a;) = /iR/mR, where /iR and mR denote the renormalized twisted and 
untwisted quark masses, w has been tuned to tt/2 by adjusting mo appropriately (cf. [32] for 
details). As argued in |24] this ensures automatic 0{a) improvement for static-light spectral 
quantities, e.g. mass differences between static-light baryons and the lightest static-light meson 
(the "S/S* meson"), the quantities we are focusing on in this work. 

The ensembles of gauge field configurations we are considering are listed in Table [TJ They 
correspond to a single value of the lattice spacing a ~ 0.079 fm, but various values of the pion 
mass in the range 340 MeV < mps <525MeV. The lattice extension is L"^ x T = 24^ x 48, which 
amounts to L ~ 1.9 fm and mpsL>3.3. Details regarding the generation of these gauge field 
configurations and computation and analysis of standard quantities (e.g. lattice spacing or pion 
mass) can be found in [32l [33] . 



p 


xT 




a in fm 


mps in MeV 


# of gauges 


3.90 


243 X 48 


0.0040 


0.079(3) 


340(13) 


200 






0.0064 




423(16) 


50 






0.0085 




485(18) 


50 






0.0100 




525(20) 


50 



Table 1: ensembles of gauge field configurations (a and mps have been taken from [33]; # of gauges: 
number of gauge field configurations considered). 
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We treat static-light baryons containing valence s quarks in a partially quenched approach, 
where the mass of these quarks, /iq, valence s = 0.022, is approximately equal to the mass of the 
physical s quark taken from a study of strange mesons using the same gauge field configurations 
[3i\ [35] • Note that partially quenched s quarks can be realized in two ways, either with a 
twisted mass term +i/Uq,vaicnce s75 or — i/iq,vaicnce s75 corresponding to the upper and the lower 
entry in the quark field doublet x respectively. We consider both possibilities and denote them 

In Table [T] we also list the number of gauge configurations, on which we have computed static- 
light baryon correlation functions. 
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3 The static-light baryon spectrum 



With static-light baryons we refer to baryons made from a single static quark and two light 
quarks, which can either be u, d and/or s. 



3.1 Static-light baryon trial states 

3.1.1 Static-light baryon creation operators in the continuum 

We start by discussing symmetries and quantum numbers of static-light baryons and correspond- 
ing creation operators in the continuum. 

The continuum analogs of our lattice static-light baryon creation operators are 

0Sa);U(x) = 6'^'"=Q'^(x)((^'''(^)(x))^CrV'^'(2)(x)), (4) 

where Q is a static quark operator and are light quark operators (in the usual physical 
basis). The upper indices a, b and c are color indices, C = 7072 is the charge conjugation matrix 
and r is a combination of 7 matrices, i.e. a 4 x 4 matrix acting in spin space. 

Since there are no interactions involving the static quark spin, it is appropriate to label static- 
light baryons by the angular momentum of their light degrees of freedom j. For creation opera- 
tors (j4]) it is determined by V and can either be j = or j = 1. j = states correspond to total 
angular momentum J = 1/2, while j = 1 states correspond to degenerate pairs of states with 
total angular momentum J =1/2 and J = 3/2, respectively. 

Parity is also a quantum number depending on T. Either V = + ox V = —. 

The flavor quantum numbers are isospin / and strangeness 5. To access all possible combi- 
nations, we consider light quark flavors ^(^^^/^^^^ = ud — du (corresponding to / = 0, 5 = 0), 
^(1)^(2) ^ S^uu , dd , ud + du} (corresponding to / = 1, 5 = 0), il^^^^ ip^'^'^ G {us, ds} (correspond- 
ing to / = 1/2, S = —1) and ip^^^ip^'^^ = ss (corresponding to / = 0, S" = —2). 

Creation operators Op^J(i)^^(2) ^^"^ quantum numbers of their associated trial states O^'^J^Jj^^^j) 1^) 
are collected in Table [2j Note that certain T , tjj^^^ ij:^"^^ combinations do not need to be consid- 
ered, since the corresponding creation operators are identical zero due to the anticommutation 
property of quark operators. Such F, ^/^(^)^(^) combinations are either omitted from the table or 
marked with "X". 

3.1.2 Static-light baryon creation operators in twisted mass lattice QCD 

Twisted basis lattice static-light baryon creation operators are of similar form, 

0™(2)(x) = 6'^^'=Q'^(x)((/'(i)(x))^CFx'^-(2)(x)), (5) 

where physical basis quark operators have been replaced by their twisted basis lattice counter- 
parts. 
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r 


■V 

] 


J 




s 


name 


/ 


5 


name 


I 


s 


name 


75 


n+ 




n 


n 




1 12 


_1 




X 


X 


X 


7o75 


0+ 


1/2 










1/2 


-1 


^b 


X 


X 


X 


1 


0~ 


1/2 










1/2 


-1 




X 


X 


X 


70 


0~ 


1/2 


1 







1/2 


-1 







-2 




7j 


1+ 


1/2, 3/2 


1 







1/2 


-1 







-2 


^b 


7o7i 


1+ 


1/2, 3/2 


1 







1/2 


-1 







-2 




7i75 




1/2, 3/2 










1/2 


-1 




X 


X 


X 


7o7i75 


\- 


1/2, 3/2 


1 







1/2 


-1 







-2 





Table 2: continuum static-light baryon creation operators and their quantum numbers (j^: angular 
momentum of the light degrees of freedom and parity; J: total angular momentum; /: isospin; 
S: strangeness; name: name of the corresponding b baryon(s) in |47j): operators marked with 
"X" are identically zero, i.e. do not exist. 



In the continuum the relation between the physical and the twisted basis is given by the twist 
rotation ■0 = exp^i^^T^u /2)x, where u = tt/2 at maximal twist. At finite lattice spacing, 
however, issues are more complicated: the twist rotation only holds for renormalized operators 
and the QCD symmetries isospin and parity are explicitely broken by 0{a). Nevertheless, it is 
possible to unambiguously interpret states obtained from correlation functions of twisted basis 
operators in terms of QCD quantum numbers as we will explain and demonstrate below. 

On the lattice rotational symmetry is reduced to symmetry with respect to cubic rotations. 
There are only five different representations of the cubic group Oh corresponding to integer 
angular momentum j. j = in the continuum corresponds to the Ai representation on the lattice 
containing angular momenta j = 0,4,7,..., while j = I corresponds to the Ti representation 
containing j = 1, 3, 4, . . . 

While in twisted mass lattice QCD the z-component of isospin 1^ is still a quantum number, 
isospin / and parity V are explicitely broken by the Wilson term, which is proportional to the 
lattice spacing. Only a specific combination of both symmetries, light flavor exchange combined 
with parity, is still a symmetry in twisted mass lattice QCD. We denote this symmetry by 
acting on the light twisted basis quark doublet x = (x^^^X^*^^) according to 'P^^'^^x = 7o''"iX) 
where ti is the first Pauli matrix acting in flavor space. Consequently, the four QCD sectors 
labeled by / = 0, 1 and V = +,— are pairwise combined. "P^*'^) = -|- is a combination of 
(7 = 0,7'=-) and {I = 1,V = +), while P^*"') = - is a combination of (/ = 0,P = +) and 
{I=1,V = -). 

As explained in section [2] the partially quenched s quark can be realized in two ways denoted by 
X^*^^ and x^^~K respectively. As a consequence baryons computed at finite lattice spacing on 
the one hand with quarks and on the other hand with s~ quarks, but which are otherwise 
identical, may differ in mass. Due to automatic 0{a) improvement of twisted mass lattice QCD 
this mass splitting, however, will only be ©(a^), i.e. is expected to be rather small and will 
vanish quadratically, when approaching the continuum limit. 

Since 7''-*™-* and Iz do not commute, they cannot simultaneously be chosen as quantum numbers. 
An exception are states with 7^ = 0, which can also be classified with respect to V^*^\ 
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The lattice static-light baryon creation operators we have been using are collected in Table [3l 
Tabledland TableO Creation operators are sorted according to the twisted mass lattice quantum 
numbers of their associated trial states, i.e. creation operators exciting states from different 
sectors are separated by horizontal lines. To interpret these twisted basis creation operators in 
terms of QCD quantum numbers, we have performed an approximate rotation to the physical 
basis (neglecting renormalization and using w = '7r/2). The resulting so-called pseudo physical 
basis creation operators together with their corresponding QCD quantum numbers are also listed 
in the tables. 



twisted basis lattice operator 


pseudo physical basis operator 


r 


^(1)^(2) 


Iz 


•p(tm) 


r 


^(1)^(2) 


I 


h 


V 


name 






Ai representation = 


J = 0,4,7,... 








75 


ud — du 







75 


ud — du 








+ 


Ab 


70 


ud + du 







7o75 


ud — du 








+ 


A;, 


7075 


ud — du 







70 


ud + du 


1 









1 


ud — du 





+ 


1 


ud — du 












70 


uu/dd 


+1/-1 


XXX 


70 


uu/dd 


1 


+1/-1 










Ti representation = 


J = 1,3,4,... 








7j75 


ud — du 





+ 


7j 


ud + du 


1 





+ 




7o7i 


ud + du 





+ 


7o7i 


ud + du 


1 





+ 




7j 


ud + du 





+ 


7i75 


ud — du 












7o7j75 


ud + du 







7o7j75 


ud + du 


1 









7j 


uu/dd 


+1/-1 


XXX 


7j 


uu/dd 


1 


+1/-1 


+ 




7o7j75 


uu/dd 


+1/-1 


XXX 


7o7i 


uu/dd 


1 


+1/-1 


+ 




7o7i 


uu/dd 


+1/-1 


XXX 


7o7j75 


uu/dd 


1 


+1/-1 







Table 3: S = lattice static-light baryon creation operators and their quantum numbers; (j: angular 
momentum of the light degrees of freedom; /: isospin; z-eomponent of isospin; V: parity; 
•p(t'n): twisted mass parity ["xxx" indicates that is not a quantum number for the 

corresponding trial state]; name: name of the corresponding b baryon(s) in |47l ). 



3.1.3 Smearing of gauge links and quark fields 

To enhance the overlap of the trial states ^r™^w^(2) *° lyiiig static-light baryon states, we 
make extensive use of standard smearing techniques. This allows to read off static-light baryon 
masses from correlation functions at rather small temporal separation, where the signal-to-noise 
ratio is favorable. 

Smearing is done in two steps. At first we replace all spatial gauge links by APE smeared 
versions. The parameters we have chosen are A^ape = 40 and oape = 0.5. Then we use Gaussian 
smearing on the light quark fields x^"\ X^'^K X^**""^^ ^ind ^ which resorts to the APE smeared 
spatial links. We consider three different smearing levels, characterized by A'^Qauss € {10 , 40 , 90} 
and KGauss = 0.5. This amounts to light quark field operators with approximate widths of 



7 



twisted basis lattice operator 




pseudo 


physical basis operator 




r 


y(l)y(2) 

A, A, 


Iz 


r 


^(1)^(2) 


/ 




V 


name 






A\ representation 


= J = 0,4, 7, 










L 


Cio / ilo 


+1/2/ - 1/2 


75 


us/ds 


1/2 


+ 1/2/ - 


1/2 


1 




/u /o 


us^ /ds~ 


+1/2/ - 1/2 


7075 


us 1 ds 


1/2 


+1/2/ - 


1/2 


+ 




1 o 


us^ /ds~ 


+1/2/ - 1/2 


1 


us/ds 


1/2 


+1/2/ - 


1/2 






Tn 


us'^ / ds~ 


+1/2/ - 1/2 


70 


us/ds 


1/2 


J-1 /9 / 
+ i/^/ - 


1 /9 






75 


lie / /V t> 

Lib 1 LLo 


+1/2/ - 1/2 


75 


us/ds 


1/2 


+1/2/ - 


1/2 


1 








+1/2/ - 1/2 


7075 


us/ds 


1/2 


+1/2/ - 


1/2 


+ 




1 




+1/2/ - 1/2 


1 


us / ds 


1/2 


+1/2/ - 


1/2 








us~ /ds^ 


+1/2/ - 1/2 


70 


us/ds 


1/2 


-1-1/9/- 
-|-i/^/ 


1 /9 










T\ representation 


= J = 1,3,4,. 










7i 


us~^/ds~ 


+1/2/ - 1/2 


7i 


us/ds 


1/2 


+1/2/ - 


1/2 


+ 




7o7i75 


us^/ds~ 


+1/2/ - 1/2 


7o7i 


us / ds 


1/2 


+1/2/ - 


1/2 


+ 




7i75 


us~^/ds~ 


+1/2/ - 1/2 


7i75 


us/ds 


1/2 


+1/2/ - 


1/2 






7o7i 


us^ /ds~ 


+1/2/ - 1/2 


7o7i75 


us/ds 


1/2 


+1/2/ - 


1/2 






7i75 


us~ /ds^ 


+1/2/ - 1/2 


I'j 


us / ds 


1/2 


+1/2/ - 


1/2 


+ 




7o7i 


us~ /ds'^ 


+1/2/ - 1/2 


7o7i 


us/ds 


1/2 


+1/2/ - 


1/2 


+ 




7i 


us~ /ds^ 


+1/2/ - 1/2 


7i75 


us/ds 


1/2 


+1/2/ - 


1/2 






707j75 


us~ /ds^ 


+1/2/ - 1/2 


707i75 


us/ds 


1/2 


+1/2/ - 


1/2 







Table 4: 5* = — 1 lattice static- light baryon creation operators and their quantum numbers; (j: angular 
momentum of the light degrees of freedom; /: isospin; 1^: z-component of isospin; V: parity; 
name: name of the corresponding h baryon in [47]). 



{1.58 X a, 3.16 x a, 4.74 x a} « {0.12 fm, 0.25 fm, 0.37 fm} (cf. [24] for details). 
Smeared static light baryon creation operators are denoted by S^'^'^'^^^iO^^^^^y^^^^^. 

3.2 Correlation matrices 

For each sector characterized by strangeness 5, angular momentum of the light degrees of free- 
dom J, z-component of isospin J^, and in certain cases twisted mass parity we compute 
temporal correlation matrices 

'^{r,,{x(i)x(2'),,A^Gauss,i),(rfc,(x(i)x(2))fc,AfGausB,fc)(*) = 

We consider all the creation operators listed in Table El Table S] and Table [5] at three different 
smearing levels -/Vcauss 

G {10, 40, 90} as explained in the previous subsection. This amounts 
dependent on the sector to 3 x 3, 9 x 9 or 12 x 12 correlation matrices. 

Static quarks are treated with the HYP2 static action [SG] [371 [38] , Wilson lines appearing in 
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twisted basis lattice operator 


pseudo physical basis operator 


r 


^(1)^(2) 


h 


r 


^(1)^(2) 


/ 


h 


V 


name 




Ai representation = 


J = 0,4,7,.. 








70 


s 





70 
















Ti representation = 


J = 1,3,4,... 








7j 







7i 


ss 








+ 




7o7j75 


s'^s'^ / s~ s~ 





7o7j 


ss 








+ 




7o7i 







7o7i75 


ss 













Table 5: S = —2 lattice static-light baryon creation operators and their quantum numbers; (j: angular 
momentum of the light degrees of freedom; I: isospin; 1^: z-component of isospin; "P: parity; 
name: name of the corresponding b baryon in j47j ) . 



static quark propagators are formed by products of HYP2 smeared temporal links (cf. [21] for 
details) . 

Light quark propagators are estimated by means of Z2 x Z2 stochastic timeslice sources (cf. 
|24j for details). On each gauge field configuration we invert 48 independently chosen sources, 
all located on the same timeslice, 12 for each of the four possible light quark propagators u, d, 
s"*" and s~ . Multiple inversions of the same timeslice of the same gauge field configuration are 
beneficial with respect to statistical precision, because each correlation function contains two 
light quark propagators. This allows to form 12 x 12 = 144 statistical samples, i.e. the number 
of samples is the square of the number of inversions (cf. [6]). 



3.3 Determination of static-light baryon masses 

From correlation matrices ^ we compute effective mass plateaus by solving generalized eigen- 
value problems 

Cit)Vn{t,to) = Xnit,to)C{to)Vn{t,to) , mf{t,to) = lu ^"^^'^"^ (7) 

^n{t + t, to) 

with to = 1 (cf. e.g. |39| , l40 j ). Instead of using the full 3 x 3, 9 x 9 or 12 x 12 correlation matrices 
we have chosen "optimal submatrices" in a sense that on the one hand effective masses exhibit 
plateaus already at small temporal separations t and that on the other hand statistical errors 
on are minimized. We found that with the following choice both criteria are adequately 
fulfilled: 



• 3x3 correlation matrices: 

use 2x2 submatrices with smearing levels iVcauss S {40 , 90}; 

• 9x9 correlation matrices: 

use 3x3 submatrices with smearing levels -/Vcauss — 

90; 

• 12 X 12 correlation matrices: 

use 4x4 submatrices with smea.ring levels A^Gauss — 

90. 
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To demonstrate the quality of our lattice results, we show in Figure [T] examples of effective mass 
plateaus (at light quark mass fiq = 0.0040) corresponding to Aj, (5 = 0, / = 0, = 0+), Of, 
{S = -2, / = 0, j'^ = 1+) and its parity partner (5 = -2, / = 0, j'^ = 




2 4 6 8 10 12 2 4 6 8 10 12 

temporal separation t temporal separation t 

Figure 1: effective masses m'^ as functions of the temporal separation t at light quark mass /iq = 0.0040; 

left: Af, (5 = 0, / = 0, j"^ = 0+) from a 3 x 3 correlation matrix; right: [S = -2, I = 0, 
~ 1+) and its parity partner (5* = —2, / = 0, 2^ ~ 1^) from a 3 x 3 correlation matrix. 

We extract static-light baryon masses by fitting constants to these plateaus in regions of suf- 
ficiently large temporal separation tmin • • • imax- We found that tmin = 4 yields reasonable 
values, which are OiV) for all states investigated, tmax on the other hand hardly affects the 
resulting static-light baryon masses (on the "imax-side" of the effective mass plateau statistical 
errors are rather large and, therefore, data points have a negligible effect on the fit). The re- 
sulting fits for the examples shown in Figure [T] are indicated by dashed lines. We checked the 
stability of all our results by varying tjam by ±1- We found consistency within statistical errors. 

To assign appropriate QCD quantum numbers to the extracted static-light baryon states, we 
follow a method introduced and explained in detail in [41], section 3.1 ("Method 1: solving 
a generalized eigenvalue problem"). For the n-th state the components of the corresponding 
eigenvector Vn.j characterize the contribution of the j-th static-light baryon creation operator 
entering the correlation matrix. After transforming these operators from the twisted basis to 
the pseudo physical basis by means of the twist rotation = exp(z75r3ti;/2)x, u = 11/2 (cf. the 
right columns of Table [3l Table H] and Tabled]), one expects and and also finds that for each 
extracted state operators corresponding to only one of the two QCD sectors corresponding to 
the investigated twisted mass lattice QCD sector clearly dominate, while the contribution from 
operators from the other sector are rather small. This allows to unambiguously assign a QCD 
label to each extracted static-light baryon state. An example, the identification of Vtb {S = —2, 
/ = 0, = 1"^) and its parity partner {S = —2, 7 = 0, = 1~), is shown in Figure [2] (cf. 
also Figured] for the corresponding effective masses both having twisted mass quantum numbers 
(5 = -2,j = 1,1 = 0)). 

Since static-light baryon masses diverge in the continuum limit due to the self energy of the static 
quark, we always consider mass differences of these baryons to the lightest static-light meson 
("S/i?* meson"). In such differences the divergent self energy exactly cancels. We take the mass 
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2 4 6 8 10 12 2 4 6 8 10 12 

temporal separation t temporal separation t 



Figure 2: eigenvector components Ifn.jP as functions of the temporal separation t and their associated 
QCD quantum numbers at light quark mass = 0.0040 corresponding to the 3x3 correlation 
matrix with twisted mass quantum numbers (S ~ —2, j = 1, / = 0) (cf. also Figure[T]); left: 
ground state identified as 0,1, {S = —2, / = 0, = 1+); right: first excited state identified 
as parity partner of fib {S = —2, / = 0, = 1^). 



values of the lightest static- light mesons from [25], where they have been computed using the 
same lattice setup. The mass differences Am^*'^*(5', I, i^)a = (m(baryon : 5, 1 ^j'^) — m{B / B*))a 
(in lattice units) together with the pion masses mpsa (also in lattice units; cf. Table [Hand |33j ) 
serve as input for the extrapolation procedure to the physical u/d quark mass described in the 
next subsection. 

3.4 Extrapolation to the physical u/d quark mass 

The mass differences Am^^''^^{S,I,j''^)a = (m(baryon : S,I,j'P) - m{B/B*))a obtained for the 
four ensembles listed in Table [H which only differ in the value of the u/d quark mass (both sea 
and valence), are plotted against (mpsa)^ in Figure [IKS' = 0, / = 0, i.e. A baryons). Figure H] 
(5 = 0, I = 1, i.e. S baryons), Figure [5] [S = —1, i.e. H baryons) and Figure [6] [S = —2, i.e. O 
baryons) and are collected in appendix 

For the extrapolation to the physical u/d quark mass one could use an effective field theory ap- 
proach (Chiral HQET for example) as used e.g. to study static-light meson decay constants [42j . 
However, this approach has not fully been developed to discuss mass differences Am'^*^*(S', /, j'^)a 
between excited static-light baryon states and the lightest static-light meson so is not appropri- 
ate here. Instead we use the simplest assumption, which is supported by our results: a linear 
dependence in (mpsa)^. 

Data points ((mpsa)^ , /^rrf^''\S,I,f)a) are correlated via (mpsa) in case they correspond to 
the same ensemble, i.e. to the same value of the u/d quark mass. We take that into account via 
a covariance matrix, which we estimate by resampling mpsa and all extracted static-light mass 
differences Am^^^^{S,I,j^)a (10000 000 samples). Consequently, we do not fit straight lines to 
the data points ((mpsa)^ , Am^*^*(5, /, j^)a) individually for every static-light baryon state, but 
perform a single correlated fit of 23 straight lines to the 23 mass differences considered. During 
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Figure 3: mass differences of S* = 0, / = static-light baryons (A baryons) to the lightest static-light 
meson as functions of (mpga)^; straight lines represent linear extrapolations to the physical 
u/d quark mass. 
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Figure 4: mass differences of S* = 0, / = 1 static-light baryons (E baryons) to the lightest static-light 
meson as functions of (mpga)^; straight lines represent linear extrapolations to the physical 
u/d quark mass; plots in the same line only differ in (left: 1^ = 0; right: 1^ — ±1). 
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Figure 5: 



mass differences of 5* = — 1 static-light baryons (S baryons) to the lightest static- light meson 
as functions of (mpso)^; straight lines represent linear extrapolations to the physical u/d 
quark mass; plots in the same line only differ in the sign of the twisted mass term of the s 
valence quark (left: us~^/ds~; right: us~ /ds~^). 
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Figure 6: mass differences oi S = —2 static-light baryons (fi baryons) to the lightest static-light meson 
as functions of (mpsa)^; straight lines represent linear extrapolations to the physical u/d 
quark mass; plots in the same line only differ in the signs of the twisted mass terms of the s 
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the fitting we take statistical errors both along the horizontal axis (errors in mpsa) and along 
the vertical axis (errors in Am^^^^{S, I,j^)a) into account. The method for performing such 
two-dimensional fits is explained in detail in 



We find that a fit, which is linear in the light quark mass (represented by the mass squared of 
the light-light pseudoscalar meson (mpga)^) is acceptable, i.e. yields x^/dof ~ 0.59 <1. This 
fit enables us to extrapolate to the physical u/d quark mass, in this work taken as mps = 
135 MeV and converted to lattice units by using the lattice spacing a = 0.079 fm [33] resulting 
in (mpgo)^ = 0.054^ (cf. Figure O Figure [H Figure O and Figure ED . 

Extrapolations of static-light mass differences to the physical u/d quark mass are listed in TableE] 
in MeV. Since there seems to be a controversy of around 10% regarding the value of the lattice 
spacing in physical units, when using on the one hand the pion mass 771,^ and the pion decay 
constant /tt and on the other hand the pion mass rrij^ and the nucleon mass tun to set the scale 
(a = 0.079(3) fm |33j versus a = 0.089(5) fm |43]). we also list dimensionless ratios of static-light 
mass differences, 

= ^T'l^'lC'" . (8) 

These ratios are pretty independent of the lattice spacing and, therefore, preferable, when making 
predictions or when comparing to other lattice or model computations or to experimental results. 

For static-light baryons with S = —1 and S = —2 our results depend on the bare s quark 
mass chosen. We use ^q,vaicncc s = 0.022 taken from studies of strange-light mesons [Ml [35] as 
mentioned in section [2j Possible systematic errors arising from a slightly incorrect value of the 
s quark mass are expected to be smaller than the corresponding statistical errors, because the 
mass differences we compute turn out to be rather weakly dependent on the masses of their light 
valence quarks (cf. Figure [3l Figure [Hand Figure [5]). A possibility to estimate their magnitude 
for is to estimate the slope of m{Qb) as a function of /Uq^ valence s by means of the experimental 
results on m(r2;,) and m{T,i,) and assuming a systematic error of 10% for the lattice spacing (i.e. 
roughly the difference obtained with the two scale setting methods j33|. |43] ) . Then one arrives 
at a systematic error of around 0.1 x (m(r2b) — m(Sb)) ~ 23 MeV. This number is consistent 
with an even simpler method of estimation namely just taking a systematic error of 0.1 x ijig 
for every s valence quark, where rUg = 80 MeV . . . 130 MeV [17|- We intend to investigate the 
s quark dependence in more detail and to quantify the corresponding systematic error more 
precisely in a subsequent publication. 

Static-light baryon states with (S* = 0, / = 1) and either Iz = {ud + du) or 1^ = ±1 {uu/dd) 
are not degenerate in twisted mass lattice QCD, but differ by discretization errors. These 
discretization errors are, however, only 0{a^) and, therefore, expected to be rather small. As 
can be seen from Table [6] = and Iz = ±1 states agree within statistical errors, which is a 
strong indication that discretization errors are indeed negligible. For the ratios R^^^^{S,I,j'^) 
and for interpolations to the physical b quark mass (cf. section H]) we subsequently use Iz = 
results. 

For (S = —1, / = 1/2) static-light baryon states similar statements apply. 

For static-light baryon states with two s quarks, i.e. {S = —2, 7 = 0), the situation is somewhat 
different. On the operator level one should not use different lattice discretizations of the two s 
quarks, i.e. one twisted and one twisted s~ quark (for example the operator with F = 75, 
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s 


/ 
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name 


flavor 


a from [33] 


a from |43| 










0+ 




ud — du 


430(40) 


382(39) 


0.480(42) 










ud — du 


1499(156) 


1330(149) 


1.672(166) 






1" 




ud — du 


681(116) 


605(106) 


0.760(127) 





1 


1+ 




ud + du 
uu/dd 


611(45) 
635(44) 


543(46) 
563(45) 


0.682(45) 










ud + du 
uu/dd 


905(176) 
986(137) 


803(160) 
876(127) 


1.010(194) 






1" 




ud + du 
uu/dd 


1037(88) 
957(117) 


921(87) 
850(110) 


1.158(91) 


-1 


1/2 


0+ 




us~ /ds^ 
us'^ /ds^ 


602(33) 
629(41) 


534(37) 
558(43) 


0.672(30) 






1+ 




us~ /ds^ 
us'^ /ds^ 


767(39) 
778(38) 


681(44) 
690(44) 


0.856(33) 










us~ /ds^ 
us~^ /ds~ 


1205(117) 
954(141) 


1070(113) 
847(130) 


1.351(123) 






1" 




us~ /ds^ 
us'^ /ds^ 


1062(71) 
1068(91) 


943(74) 
948(89) 


1.185(69) 


-2 





1+ 




s'^s'^ /s~s~ 
s~^s^ /s~s^ 


896(39) 
896(38) 


795(48) 
795(47) 


1 










s^s'^ /s^s^ 
s'^s^ /s~s^ 


1336(64) 
1296(94) 


1186(75) 
1150(96) 


1.491(53) 






1" 




s~^s~^ /s~s~ 
s~^s~ /s~s'^ 


1236(76) 
1255(61) 


1097(81) 
1114(71) 


1.380(72) 



Table 6: static-light mass differences Am"*^'(S', /, j^) = TO(baryon : S,1 ,3'^) - m{B / B*) in MeV (scale 
setting via m^r and a = 0.079(3) fm [33], and via m^r and toat, a = 0.089(5) fm [43]) and 
dimensionless ratios of static-light mass differences (cf. ([8])). 
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t^(i)t^(2) = gg — gg clearly is identically zero, while T = 75, x^'^^X^'^^ = s~^s~ — s~s~^ would 
give a non-zero correlator, because there is no propagation from to s^). Only on the level 
of correlators one can insert one s"*" lattice propagator and one s~ lattice propagator in a 
meaningful way. The corresponding masses are also listed in Table [6] and agree with their 
s~^s~^/s~s~ counterparts. For the ratios W^^^{S, I , j^) and for interpolations to the physical h 
quark mass (cf. section we subsequently use s"^s+/s^s" results. 

In principle contamination of some of the computed static-light baryon masses by multiparticle 
states (e.g. by a lighter static-light baryon and a pion) can at this stage not rigorously be 
excluded. However, from previous experience with similar lattice computations one strongly 
expects that the overlap of the used trial states to multiparticle states is extremely small and, 
therefore, that contamination of the obtained static-light baryon masses is negligible within 
statistical errors. A possibility to investigate this issue in detail is to compute matrix elements 
of two operators, where one is a "single particle baryon operator" and the other a "multiparticle 
operator" (cf. [Ml HU |46] , where the method has been introduced and applied to glueballs and 
various types of mesons). Using this method we were e.g. able to confirm that the parity partner 
of the lightest static- light meson is essentially not affected by multiparticle states (cf. |25|). 

Finally it is interesting to compare our static-light mass differences to recent results obtained by 
other lattice groups. In [8] three h baryon states. A;,, S^/S^ and O;,, are computed. The method 
of scale setting used in this paper, imposing ro = 0.49 fm, can easily be applied to our lattice 
results, which are then in excellent agreement within statistical errors. Since in and \lQi\ ITT] 
scale setting methods are used, which are less straightforward to adapt (via mir, rriK and 
and chiral perturbation theory [9] and via T lllj). we directly compare the resulting mass 
differences in MeV. When comparing to those of our results corresponding to a = 0.089(5) fm 
(i.e. scale setting via mTv), we also find agreement for all states computed in [9l [ini [II], A;,, 
Sb/S*, Hft, H', {S = -2, I = 1/2, f = 1+) and 0,,. 

We are also able to predict eight static-light baryon states of negative parity, for which no other 
lattice results seem to be available at the moment. Therefore, we compare these predictions to 
the quark model calculation in [23] . Also here we find agreement within statistical errors. 
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4 Interpolation to the physical b quark mass 



To make contact with experimentally available results on the spectrum of b baryons, we need 
to correct for the finite mass of the b quark. In Heavy Quark Effective Theory the leading 
correction is ©(l/mg), where uiq is the mass of the heavy quark. It is possible in principle to 
compute such corrections from first principles by means of lattice QCD (cf. e.g. j2H 122]). This 
we intend to explore in the future, but here we use a more direct method, to establish the size 
of the correction between static quarks and b quarks of physically realistic mass. 

We linearly interpolate in mc/mq between our static-light lattice results and corresponding 
experimental data for charmed baryons. As a measure of the heavy quark mass mq we take 
the masses of the ground state heavy-light mesons {D or S), i.e. we interpolate to rric/mi, = 
m{D)/m{B) = 0.35. This measure is equivalent to another (such as using quark masses in some 
scheme) to the order l/mg we are considering. One test of this interpolation can be made: 
the hyperfine splitting between Sc(2520) and Sc(2455) is around 64MeV; interpolating between 
this number and the static limit results in 0.35 x 64MeV = 22MeV, which is in fair agreement 
with the observed splitting of m(S^) — m(Sb) = 21MeV [47]. Results of these interpolations 
are collected in Table [71 Note that these rric/raQ corrections break the heavy spin degeneracy 
of static- light baryon states with j = 1. 

4.1 Discussion of possible systematic errors 

Our lattice results might be associated with certain systematic errors, which we list and briefly 
discuss in the following. 

• Scale setting: 

the dominating source of systematic error arises from the ambiguity introduced by the 
two methods of scale setting \33\ 143] . which is around 10%. Although it seems that the 
lattice spacing a = 0.089(5) fm determined by means of the nucleon mass seems more 
appropriate, when comparing to experimental results or to publications from other lattice 
collaborations (cf. section 13.41 and section 14. 2p . we strongly recommend to consider the 
dimensionless ratios i?^^*'*, where scale setting errors are essentially eliminated. 

• Extrapolation to the physical u/d quark mass: 

as explained in section[2]our results, which cover pion masses in the range 340 MeV < mps ^ 525 MeV, 
are consistent with a linear dependence in (mps)^ for all static- light baryon states. Whether 
there are deviations at significantly lighter u/d quark masses, will be studied using corre- 
sponding ETMC gauge field configurations, which will be available soon. 

• Possibly incorrect tuning of the s quark mass: 

this issue has already been discussed in section 13.3] where we estimate the systematic 
error for static-light baryons with a single valence s quark to be around 10 MeV and for 
those with two valence s quarks to be around 20 MeV. Note that the extrapolation to the 
physical b quark mass by means of experimental results on c baryons reduces these errors 
by around 1/3. 

• Extrapolation to the physical b quark mass: 

the validity of the interpolation between static lattice results and charm experimental 
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Table 7: 6 baryon mass differences Am{S, I, J^) = TO(baryon : S*, /, J^) — m{B) in MeV (scale setting 
via m^r and /jr, a = 0.079(3) fm |33| . and via and m^r, a = 0.089(5) fm [43]) and dimen- 
sionlcss ratios of baryon mass differences differences (cf. (jS])); ^ experimental results have been 
taken from |47j with exception of m(51b), which has been taken from [5]; ^ this number docs not 
require any lattice result. Lines associated with quantum numbers, where no corresponding c 
baryons have experimentally been measured are filled with — . 
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results by the order l/mq of HQET has been tested for baryons via S^/S^ as explained 
above and for mesons via B/B* [25]. These two tests indicate validity up to ~ 5%. Since 
the hyperfine splitting in the b region is of order 20MeV to 50MeV (cf. e.g. B/B*, B1/B2, 
Y^b/'^l in [13) one expects a corresponding systematic error of <2.5MeV. 

• Electromagnetic and isospin breaking effects: 

experimental results on and indicate that such effects can be of order 5MeV to 
lOMeV. 

• Neglect of s and c sea quarks: 

the systematic error arising from our neglect of the s and c quark contribution to the sea 
is expected to be significantly smaller than current statistical errors. Will will address 
and quantify this error in a future study making use of recently generated Nj = 2 + 1 + 1 
ETMC gauge field configurations [igjlig]. 

• Continuum limit: 

since we use a rather fine lattice spacing and an 0{a) improved lattice formulation (twisted 
mass lattice QCD at maximal twist), we expect discretization effects to be negligible. This 
expectation is supported by the computation and comparison of / = 1 states (/^ = ±1 
versus Iz = 0), which are not degenerate in twisted mass lattice QCD, but differ by 
O(a^). This constitutes a direct check of lattice discretization effects, for which we found 
no indication (cf. Table [6]). Moreover, we have recently investigated the continuum limit 
for b mesons using the same gauge field configurations and also did not find any sign of 
discretization effects |25j . 

• Multiparticle states: 

contamination of static-light baryon states by multiparticle states of the same quantum 
number have been discussed in section 13. 3i It seems rather unlikely that they introduce a 
systematic error, which is significant compared to current statistical errors. 

In total the sum of these systematic errors should not exceed 25 MeV, which is of the same order 
of magnitude as for our recent lattice results on B mesons [25], where we quoted a maximal 
systematic error of 20 MeV. An additional uncertainty of 10% should be assigned, when consid- 
ering mass differences in MeV, i.e. Am'^*'^(S', I, j'^). For the dimensionless ratios R^'^^'^{S,I,j'^) 
collected in Table [7] the latter is not present, while the above mentioned 25 MeV translate to a 
systematic error of around 5%. 

4.2 Comparison to experimental results 

In experiments five b baryon states have been measured: A^, Sf,; 5]^, Ef, and ilf,. We compare our 
lattice results with these experimental results in Table [71 As already mentioned in the previous 
section the lattice spacing depends to some extent on the observables used to introduce physical 
units. While setting the scale via and /^r [33] yields lattice results, which are around 10% to 
20% larger than their experimental counterparts, using and ttt-tv [13] leads to significantly 
better agreement. To reduce scale setting effects as much as possible, we prefer to compare the 
dimensionless ratios R^^^'^{S,I,j^) and R'^^^'''{S, , I), which have been defined in ([8]). While 
i?®^P'^(5, /, j'^) denotes the ratio of experimentally measured b baryons, R^^^'^{S, I, j^) is the 
linear nic/mQ interpolation between the static-light ratio from Table [6] and the corresponding 
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ratio of experimentally measured c baryons. As can be seen from Table [71 there is reasonable 
agreement between our lattice and experimental results for the four available ratios. 

There are seven more b baryon states we are able to predict, but which have not yet been 
measured by experiment. Their values in MeV as well as the ratios B}^^'^{S, I, j'^) are also 
collected in Table [71 
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5 Conclusions 



We have computed the spectrum of static-hght baryons by means of lattice QCD using Nf = 2 
flavors of light quarks. We have considered all possible combinations of two light quarks, i.e. 
A, S, H and O baryons, angular momentum/spin of the light degrees of freedom j £ {0 , 1} and 
both parity V = + and V = —. In particular we were able to predict a number of negative 
parity states, which have at the moment neither been measured experimentally nor previously 
been computed on the lattice. 

We have employed the assumption of a l/mq dependence on the heavy quark mass together 
with experimental results for c baryons to allow us to estimate the spectrum that one would 
obtain for b quarks of finite physical mass. 

The wide variety of computed states (both static-light baryons and b baryons) will be a valuable 
resource for model builders and might give input for future experiments. 

Obvious directions to continue this research include (i) investigating the continuum limit; (ii) per- 
forming similar computations at lighter u/d quark masses; (iii) studying the dependence of H 
and baryons on the s quark mass; (iv) extending these computations to Nj = 2 -|- 1 -|- 1 
flavor ETMC gauge field configurations [IHl HH] ; (v) considering non-trivial gluonic excitations 
allowing to study total angular momentum of the light degrees of freedom j > 1; (vi) replacing 
experimental input for c baryons by corresponding lattice results with heavy quarks of finite 
mass [50] and/or combining such results with a recently proposed method for lattice B physics 
[5T] to compute the spectrum of b baryons in an alternative way. 
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0+ 


us~ /ds^ 


0.2419(54) 


0.2346(121) 


0.2356(77) 


0.2444(77) 








us'^/ds" 


0.2560(76) 


0.2663(86) 


0.2628(101) 


0.2671(105) 






0- 


us~ /ds^ 


0.4559(247) 


0.4065(445) 


0.4112(537) 


0.4048(298) 








us~^ /ds~ 


0.4118(320) 


0.4139(355) 


0.4130(506) 


0.4879(114) 






1+ 


us~ /ds^ 


0.3107(53) 


0.3198(79) 


0.3120(91) 


0.3203(104) 








us^ /ds^ 


0.3131(48) 


0.3119(123) 


0.3066(99) 


0.3228(117) 






l- 


us^ /ds^ 


0.4399(122) 


0.4772(242) 


0.5113(308) 


0.4568(206) 








us^ /ds^ 


0.4554(177) 


0.4666(312) 


0.5134(349) 


0.4923(275) 


-2 







s^s~^ /s~s~ 


0.5195(90) 


0.5070(143) 


0.5198(122) 


0.4879(114) 










0.4887(176) 


0.4927(364) 


0.4790(455) 


0.4397(260) 






1+ 


s'^s^ /s^s~ 


0.3508(34) 


0.3488(86) 


0.3357(80) 


0.3422(94) 








s^s~ /s~s~^ 


0.3513(35) 


0.3488(64) 


0.3349(68) 


0.3451(80) 






1- 


s~^s~^ /s~s~ 


0.5150(117) 


0.5177(287) 


0.5650(235) 


0.5281(272) 








s'^s' /s~s'^ 


0.5165(75) 


0.5300(219) 


0.5460(138) 


0.5279(180) 



24 



Acknowledgments 

It is a pleasure to thank Vladimir Galkin and Chris Michael for many hours of helpful discus- 
sions. We acknowledge further useful discussions with Jaume Carbonell, William Detmold, Karl 
Jansen, Andreas Kronfeld and Michael Miiller-Preussker. 

The major part of computations has been performed at the PC farm at DESY Zeuthen. We 
thank DESY as well as its staff for technical advice and help. 

This work has been supported in part by the DEC Sonderforschungsbereich TR9 Computer- 
gestiitzte Theoretische Teilchenphysik. 

References 

[1] T. Aaltonen et al. [CDF Collaboration], "Eirst observation of heavy baryons and S^," 
Phys. Rev. Lett. 99, 202001 (2007) [ arXiv:0706. 3868 [hep-ex]]. 

[2] V. M. Abazov et al. [DO Collaboration], "Direct observation of the strange h baryon HjJ'," 
Phys. Rev. Lett. 99, 052001 (2007) [ arXiv:0706.1690 [hep-ex]]. 

[3] T. Aaltonen et al. [CDE Collaboration], "Observation and mass measurement of the baryon 
Phys. Rev. Lett. 99, 052002 (2007) [arXiv:0707.0589 [heo-ex]]. 

[4] V. M. Abazov et al. [DO Collaboration], "Observation of the doubly strange h baryon ^^j^," 
Phys. Rev. Lett. 101, 232002 (2008) [arXiv:0808.4142, [heo-ex]]. 

[5] T. Aaltonen et al. [CDE Collaboration], "Observation of the Jl^-baryon and measurement of 
the properties of the H;,- and Ob-baryons," Phys. Rev. D 80, 072003 (2009) [arXiv:0905.3123] 
[hep-ex]]. 

[6] C. Michael and J. Peisa [UKQCD Collaboration], "Maximal variance reduction for stochas- 
tic propagators with applications to the static quark spectrum," Phys. Rev. D 58, 034506 
(1998) [arXiv:hep-lat/9802015| . 

[7] W. Detmold, K. Orginos and M. J. Savage, "Si? potentials in quenched lattice QCD," 
Phys. Rev. D 76, 114503 (2007) |arXiv:hep-lat/0703009] . 

[8] T. Burch, C. Hagen, C. B. Lang, M. Limmer and A. Schafer, "Excitations of single-beauty 
hadrons," Phys. Rev. D 79, 014504 (2009) [arXiv:0809.1103l [heo-lat]]. 

[9] W. Detmold, C. J. Lin and M. Wingate, "Bottom hadron mass splittings in the static limit 
from 2+1 flavour lattice QCD," Nucl. Phys. B 818, 17 (2009) |arXiv:0812. 25831 [hep-lat]]. 

[10] H. W. Lin, S. D. Cohen, N. Mathur and K. Orginos, "Bottom-hadron mass splittings from 
static-quark action on 2-M-flavor lattices," Phys. Rev. D 80, 054027 (2009) [arXiv:0905.4T20l 
[hep-lat]]. 

[11] H. W. Lin, S. D. Cohen, L. Liu, N. Mathur, K. Orginos and A. Walker-Loud, "Heavy- 
Baryon Spectroscopy from Lattice QCD," Comput. Phys. Commun. 182, 24 (2011) 
[arXiv:1002.4710, [hep-lat]] . 



25 



[12] M. Neubert, "Heavy quark symmetry," Phys. Rept. 245, 259 (1994) 
|arXiv:hep-ph/9306320l . 

[13] T. Mannel, "Heavy-quark effective field theory," Rept. Prog. Phys. 60, 1113 (1997). 

[14] R. Lewis and R. M. Woloshyn, "Bottom baryons from a dynamical lattice QCD simulation," 
Phys. Rev. D 79, 014502 (2009) |arXiv:0806.4"783l [heo-lat]]. 

[15] H. Na and S. Gottlieb, "Heavy baryon mass spectrum from lattice QCD with 2+1 dynamical 
sea quark flavors," PoS LATTICE2008, 119 (2008) |arXiv:0812.1235] [hep-lat]]. 

[16] S. Meinel, W. Detmold, C. J. Lin and M. Wingate, "Bottom hadrons from lattice QCD 
with domain wall and NRQCD fermions," PoS LAT2009, 105 (2009) [arXiv:0909.3837] 
[hep-lat]]. 

[17] B. A. Thacker and G. P. Lepage, "Heavy quark bound states in lattice QCD," Phys. Rev. 
D 43, 196 (1991). 

[18] R. Lewis, "Bottom and charmed hadron spectroscopy from lattice QCD," larXiv: 1010.08891 
[hep-lat] . 

[19] M. Bochicchio, G. MartinelU, C. R. Allton, C. T. Sachrajda and D. B. Carpenter, "Heavy 
quark spectroscopy on the lattice," Nucl. Phys. B 372, 403 (1992). 

[20] D. Guazzini, H. B. Meyer and R. Sommer [ALPHA Collaboration], "Non-perturbative 
renormalization of the chromo-magnetic operator in heavy quark effective theory and the 
B*-B mass splitting," JHEP 0710, 081 (2007) [ arXiv:0705.1809. [hep-lat]]. 

[21] B. Blossier, M. della Morte, N. Garron and R. Sommer, "HQET at order 1/m: I. 
Non-perturbative parameters in the quenched approximation," JHEP 1006, 002 (2010) 
jar Xiv: 1001 .4783 [hep-lat]]. 

[22] B. Blossier, M. Della Morte, N. Garron, G. von Hippel, T. Mendes, H. Simma and R. Som- 
mer [Alpha Collaboration], "HQET at order 1/m: II. Spectroscopy in the quenched ap- 
proximation," JHEP 1005, 074 (2010) |arXiv:1004.266n [hep-lat]]. 

[23] D. Ebert, R. N. Faustov and V. O. Galkin, "Masses of excited heavy baryons in the rela- 
tivistic quark model," Phys. Lett. B 659, 612 (2008) |arXiv:0705.2957] [hep-ph]]. 

[24] K. Jansen, C. Michael, A. Shindler and M. Wagner [ETM Collaboration], "The static-light 
meson spectrum from twisted mass lattice QCD," JHEP 0812, 058 (2008) [arXiv:0810.l843] 
[hep-lat]]. 

[25] C. Michael, A. Shindler and M. Wagner [ETM Collaboration], "The continuum limit of the 
static-light meson spectrum," JHEP 1008, 009 (2010) [arXiv: 1004.4235) [hep-lat]]. 

[26] M. Wagner and C. Wiese [ETM Collaboration], "The Spectrum of static-light baryons in 
twisted mass lattice QCD," PoS LATTICE2010, 130 (2010) jarXiv:1008.0653l [hep-lat]]. 

[27] P. Weisz, "Continuum limit improved lattice action for pure Yang-Mills theory. 1," Nucl. 
Phys. B 212, 1 (1983). 



26 



[28] R. Frezzotti, P. A. Grassi, S. Sint and P. Weisz [Alpha collaboration], "Lattice QCD with 
a chirally twisted mass term," JHEP 0108, 058 (2001) [ arXiv:hep-lat/0101001] . 

[29] R. Frezzotti and G. C. Rossi, "Chirally improving Wilson fermions. 1. 0{a) improvement," 
JHEP 0408, 007 (2004) |arXiv:hep-lat /0306014| . 

[30] R. Frezzotti and G. C. Rossi, "Chirally improving Wilson fermions. 2. Four-quark opera- 
tors," JHEP 0410, 070 (2004) |arXiv:hep-lat/04 07002] . 

[31] A. Shindler, "Twisted mass lattice QCD," Phys. Rept. 461, 37 (2008) |arXiv:0707.4093] 
[hep-lat]]. 

[32] P. Boucaud et al. [ETM Collaboration], "Dynamical twisted mass fermions with light 
quarks: simulation and analysis details," Comput. Phys. Commun. 179, 695 (2008) 
[ arXiv:0803.0224 [hep-lat]]. 

[33] R. Baron et al. [ETM Collaboration], "Light meson physics from maximally twisted mass 
lattice QCD," JHEP 1008, 097 (2010) [ arXiv:0911.5061 [heo-lat]]. 

[34] B. Blossier et al. [ETM Collaboration], "Light quark masses and pseudoscalar decay con- 
stants from Nf = 2 Lattice QCD with twisted mass fermions," JHEP 0804, 020 (2008) 
|arXiv:0709.4574] [hep-lat]] . 

[35] B. Blossier et al. [ETM Collaboration], "Pseudoscalar decay constants of kaon and D- 
mesons from Nf = 2 twisted mass lattice QCD," JHEP 0907, 043 (2009) [arXiv: 0904^0954] 
[hep-lat]]. 

[36] A. Hasenfratz and F. Knechtli, "flavour symmetry and the static potential with hypercubic 
blocking," Phys. Rev. D 64, 034504 (2001) |arXiv:hep-lat /0103029] . 

[37] M. Delia Morte et al, "Lattice HQET with exponentially improved statistical precision," 
Phys. Lett. B581, 93, (2004) |arXiv:hep-lat/0307021 ]. 

[38] M. Delia Morte, A. Shindler and R. Sommer, "On lattice actions for static quarks," JHEP 
0508, 051 (2005) |arXiv:hep-lat/0506008] . 

[39] C. Michael, "Adjoint sources in lattice gauge theory," Nucl. Phys. B 259, 58 (1985). 

[40] B. Blossier, M. Delia Morte, G. von Hippel, T. Mendes and R. Sommer, "On the generalized 
eigenvalue method for energies and matrix elements in lattice field theory," JHEP 0904, 
094 (2009) |arXiv:0902.1265] [hep-lat]]. 

[41] R. Baron et al. [ETM Collaboration], "Computing K and D meson masses with A'^^ = 2-|-l-|- 
1 twisted mass lattice QCD," Comput. Phys. Commun. 182, 299 (2011) [arXiv: 1005.2042] 
[hep-lat]]. 

[42] B. Blossier et al. [ETM Collaboration], "/^ and with maximally twisted Wilson 
fermions," PoS LAT2009, 151 (2009) |arXiv:0911.3757l [hep-lat]]. 

[43] C. Alexandrou et al. [ETM Collaboration], "Nucleon electromagnetic form factors in twisted 
mass lattice QCD," larXiv:1102.2208l [hep-lat]. 



27 



[44] C. McNeile and C. Michael [UKQCD Collaboration], "Mixing of scalar glueballs and flavour- 
singlet scalar mesons," Phys. Rev. D 63, 114503 (2001) [arXiv:hep-lat /0010019] . 

[45] C. McNeile, C. Michael and P. Pennanen [UKQCD Collaboration], "Hybrid meson decay 
from the lattice," Phys. Rev. D 65, 094505 (2002) |arXiv:hep-lat /0201006] . 

[46] C. McNeile, C. Michael and G. Thompson [UKQCD Collaboration], "Hadronic decay of a 
scalar B meson from the lattice," Phys. Rev. D 70, 054501 (2004) |arXiv:hep-lat /0404010] . 

[47] K. Nakamura et al. [Particle Data Group], J. Phys. G 37, 075021 (2010). 

[48] R. Baron et al., "Light hadrons from lattice QCD with light {u,d), strange and charm 
dynamical quarks," JHEP 1006, 111 (2010) jarXiv:1004.5284l [heo-lat]]. 

[49] R. Baron et a/., "Light hadrons from Nf = 2+1 + 1 dynamical twisted mass fermions," 
PoS LATTICE2010, 123 (2010) |arXiv:1101.05T8] [hep-lat]]. 

[50] M. Papinutto, J. Carbonell, V. Drach and C. Alexandrou, "Strange and charmed baryons 
using Nf = 2 twisted mass QCD," PoS LATTICE2010, 120 (2010) [arXiv: 1012.27861 
[hep-lat]]. 

[51] B. Blossier et al. [ETM Collaboration], "A proposal for S-physics on current lattices," 
JHEP 1004, 049 (2010) |arXiv:0909.3"T87l [hep-lat]]. 



28 



